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I. INTRODUCTION 

Conventional nethods usea to represent uncerwater sounc 
transmission aave dericted transmission loss (43) as a func- 
tion of range in either tabular or graphical forn. Since 
the interest sas in determining quantities sucn as range of 
bottom trounce, convergence zones, and probability of 
detection, these nethods proved very satisfactory. However 
when it becaue necessary tc compare various underwater 
acoustic propagation loss models with each other or actual 
in-situ experiments, they provided the anaiyst with iess 
than adeguate insight into the causes for the okserved 
errors in the various aodels outputs. Furthermwore, interest 
has keen generated in analyzing the receivei signal for 
information from the sound source which cannot easily be 
induced fron transmission loss as a function or range. 

Cne method, the Wavenanber Tecanigue (WI), described by 
Fe. Re DiNapoli of NUSC [Ref. 1] and applied Ly Richard Lauer 
of NORDA [{Bef. 2], provides more information for model 
comparisons ani may have tne capability or determiaing sound 
source depth and range from the receiver in certain cases. 
DiNapoii used an analysis of acoustic propagation in the 
wavenumber domain as an internediate step in getting the 
transmission loss curve in the fast Field Program (!F?P). 
Lauer studied the wavenuaber domain information irom the FFP 
and concluded that source localization aignt be vossitle. 
Be. 3. Stamey {[Ref. 3] investigated the potential for deter- 
Mining tke source depth by using the 3rock Split~-Step Fast 
Fourier Transform (SSFFT) [Ref. 4] parabolic eguation aodel 
with toe "Wavenumper Technigue" (aT) model. His preliminary 
investigation involved the use of isospeed ani ASTREX sound 
speed proiilas with a fuliy avosorbiny bottot, varied 
source/receiver combinations, and aultipie frequencies. 
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The present investigation uses two acoustic jarakoiic 


: equation zodels for comparison with theoretical Llovd nirror 
U depictions, and the effects of observed inconsistencies ia 
: predicted transmission loss curves on the calcuiated results 
: from application of the WT model. An analysis of the sensi- 


tivity in determining source depth based on nodel output is 
attempted to ascertain possible use in naval operatiors. 
‘ To facilitate a straight-forward analysis of the results 


* and evaluation of the WT technigue for the cases studied, 
q only a selected number 9f the various computer rurs are 
4 included. Joly those curves whicn best illustrate the 


effects revealed py this research are vresented. 
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II. WAVENUMBER TECHNIQUE (WT) 


A. BACKGROUND (LLOYD MIRROR) 


The WI can be elucidated with the help of the classical 
Lloyd mirror effect [Ref. 5] which describes the interaction 
between direct path and surface reflected sound signals fron 
the same continuous wave (CW) source. Figure 2.1 illus- 
trates the geometry involved in tke Lloyd mirror effect. 


ee 


pe 


Figure 2.1 Lloyd Mirror Geometry. 
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The ccmplex pressure (P) can te expressed by 


psal -ika ee 


(4) 


when ©, and r, are related by 


(eqn 2.1) 


hs [Re (a-hy a Re «(deh)* (egn 2.2) 


with the assumption r;~ Ir, the pressure amplitude (p) can be 


approximated by 


-2fhd m=123... 2.3 
Rm mC 1E19, (eqn ) 


The interaction between the direct path and surface 
reflected waves produces constructive and destructive inter- 
ference which is manifested as peaks and nulls for the in 
phase and out of phase conditions respectively. Figure 2.2 
was produced by using equation 2.3. The null locations can 
be oktained by the relationship. 


p= 24 sin(4he ) (egn 2.4) 


When equation 2.3 is expressed in the horizontal wavenumber 
domain (k,-), the details of the resultant functional depen- 
dency on horizontal wavenumber yield useful iaformation. 
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B. PHYSICAL PROCESS 


The WT is a process by which 
(with range dependent amplitude 
into the spectral density as a 
wavenumber. The use of the AT 
ment would reguire a quadrature 
signal to attain the complex 
demodulation process can best be 
in figure 2.3 [Ref. 6]. 


2cos (ast) 


Acos (wt+x 


-2sin (at) 


Figure 2.3 Quadrature Deso 


filter 


the complex pressure wave 
and phase) is transforaned 
Eunction of the horizontal 
in the operational environ- 
demodulation of the source 
pressure. The quadrature 
defined py the illustration 


Acos((w -G)t+ x] 


Asin[(w -a)t + x] 


dulation of a Signal. 


When acoustic models are used to simuiate the environment, 


the requirement to perform a 
elizinated since the complex 


guadrature demodulation is 
pressure is directly 


a Re ee ee ee 


accessible. In actuality, acoustic parabolic equation 


models provide the complex envelope function (J) which aust 
be multiplied by the Hankel function in order to oktain the 
complex fressure. The acoustic wave is written in the fora 


P=UHS (k r)e Jot (eqn 2.5) 


where the reference wavenumker is defined by 


G) 
Ko= Co (eqn 2.6) 
and 
U=U(r, z) (eqn 2.7) 


is the solution to the appropriate parabolic equation. The 


complex pressure is corrected for voluze attenuation and 
Fourier transformed to attain the Spectral density. A pilot 
Cisplaying the results graphically with the horizontal wave- 
number on the x-axis and the normalized spectral density on 
the y-axis is constructed and the spacings between the nulls 
measured. Figure 2.4 [from Ref. 2] is an example of the iT 
for a source and receiver in the sae type of water mass. 

Since the Lloyd sirror field in tae wavenumber domain is 
given by [Ref. 2] 


Ca 


F(k) = sin( 824) e MMBZr) (eqn 2.8) 
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an alternative formulation [Ref. 2], which produces evenly 
spaced nulls, is obtained by converting the horizontal wave- 
humber to beta 


peV/x3 -&? (eqn 2.9) 


By using beta instead of the horizontal wavenumber, the 
distance, delta beta, between any two nulls can be used to 
ascertain the source depth by 


lay 2 (eqn 2.10) 


as illustrated in figure 2.5. 

The application for implying tne method of images, in 
isospeed cases, to describe the erfects of a perfectly fiat 
pressure release boundary can be justified by inspection. 
Therefore the parabolic equation can be used to produce the 
Lloyd's mirror effect and surface jlecoupling is not an 
issue, cf. Chapter 4 of Brekhovskikh [Ref. 7]. 

When the WI is implemented, care must be exercised to 
ensure that the complex pressure contains an adequate number 
of points to describe the Shortest periodicity in the pres- 
sure field. This will preclude the possibility of aliasing 
in the wavenumber domain. The real and imaginary elements 
of the pressure signal aust be modified so that the signal 
begins and ends with zero values. This modification is 
necessary because the pressure signal must represent a 
repetitive oscillation for the transform. The complex pres- 
sure array is zero-filled beyond the data in order to 
generate an array which has the length of a power of 2 for 
the Fast Fourier Iransform (FFT). The norizontal wavenumber 
increment is generated within the computer code by 
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III. ACOUSTIC MODEL ANALYSIS 


Aw PARABOLIC EQUATICNH 


The elliptical wave equation can be approximated by a 
paraktolic equation (PE) when it is assumed that the envelope 
function varies slowly with range. A detailed mathematical 
description of the PE from an acoustic point of view is 
presented by DeSanto [Ref. 8 and 9] and a simpler but less 
general description can be found in Coppens [Ref. 10]. The 
parabolic equation has the fora 


2 
oy + 2ik( SY) + (k2-K2)U =O (ean 3-1) 


where it is assumed that the pressure has the form of equa- 


tion 3.2, 


p =U(r,z)S(r) (eqn 3.2) 


S(r) represents the primary radial dependence of the field 
in terms of an outward propagating cylindrical wave of the 
form [Ref. 4] 


( 


a\ 
S(r)=eAy (kor) (eqn 3.3) 


If it is assumed that the range of interest is many wave- 
lengths from the source; then the asymptotic form of the 
Hankel function (equation 3.4) can be used (Ref. 11]. 
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‘ ss = «at 
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. HS") . \ —o e ker -1/4) ; Kor >> 1 (egn 3.4) 


Two PE models, the Brock version of the Split-Sterp Fast 

Fourier Transform (SSFFT) and toe Jeajer version of the 
Implicit Finite Difference (IFD) were investigated for 
a compatibility with the WT in deterazination of source depth. 


Be SSFFT 


1. Transaission loss Comparison with Lloyd Mirror 


To analyze sound waves in the wavenumber domain, we 
: needed a godel to frovide the pressure as a function of 
range. Initially the SSFFT nodel [Ref. 4] was selected. 
This model generates successive values for U as a function 
& of range with the helr of the algoritha 

. U(rear,z) =e darke(n?-1)/2 


4 (eqn 3.5) 

‘ FET) | 2 

{ eJArk2ke FET U(r, 2)) 

) ; This model was a natural choice since it was the model used 
S in the preliminary study dealing with the deteraination of 


the source depth by the WT [Ref. 3]. However, it soon 


x became apparent that this model's inherent weaknesses would 
“ require the consideration of another nodel if ocean botton 
- interactions were to be studied. These weaknesses, and 
~ their impact upon the WT, will be discussed shortly. 

om The SSFFT is a range dependent acoustic wave nodel 
S- which, for this analysis, will be operated in a range inde- 
4 pendent manner. In other words, only one sound speed 
he profile will be used, the water mass will be assumed homoge- 
% neous, and the bottom flat. To facilitate the study of the 


SSFFT model, the source code of the aodel was modified to 
generate the variables required for follow-on programs to 
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transforn tne pressure information and tnen display tn 


results graphically. A copy of tae source code Listing zor 
each program is provided in appendices A, 5, acdc. Taese 
programs were then linked by Job Control Language (JCL) 30 
that the WT grapnics were an autoaatically executei. 

A test cun of the SSFFI was aade for 19) az, 590 a 
source and receiver depths, 0.019 ka ranje step, faliy 
absorbing bottom, 2000 m water depth, ana an isospeei 
profile at 1590 na/sec. These values were chosen so that a 
comparison could be made with the tn2oreticai Lioyd airros 
transmission loss (figure 2.2). Tae output was disrlayed as 
a transzission loss curve (figure 3.1). From equation 2.4, 
the lccation 9£f null number 1 (R1) snaould be 33.33 xa ani 
null number 2 (K2) should be at 16.67 xn. The SSFFTI placea 
Ril at approximately 38.3 km and k2 at 17.1 ka. Since ths 
model kept the frequency, source depta, and receiver ¢c2sta 
constant, then this equates to a refareace sound speei of 
1305 and 1462 m/sec respectively. It is apparent that as 
the range decreases the error is reducei which indicates tae 
Fossibility of a problem in tne crange step section of the 
model which leads tc the introduction of a systematic and 
increasing deviation of the calculated sound field from that 
predicted ky the classical iloyits aiarror interference 
pattern. Ta2 "washing out" of the interference pattern for 
tanges less tian about 2 km is a result of the particular 
approximation of a point source used to icitiate the 
progran, which is tc be expected. Ratner rejuiar gfluctua- 
tions dévelop beyond about 30 ka; rroa tne geonetry of the 
case, it is plausible tiat taois arises froaw interzerence 
witno bottom reflected signals. 
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Anotner test was made using the same inputs as asove 2icert 
that the source and receiver were aoveé to a deptn of 1535 a 
(trigure 3.2). Since the half beamwiith is 29 degrees ani 
the bottom is fully absorbiny, tae snortest range at whicn a 
suctace ceflection ccuii occur is approximately 8.2 ka; 
beyond this range, surrtace reflections should be spaced ia 


accordaace with equation 2.4. The nulls beyond @.2 ka do 
occur at the correct ranges; however, there are slow todula- 
tions of the overall transmission loss signal. These siow 


modulations are another indication that the botton i 
fully atksorbing as nad been specified as innut to the 


program. 
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2. 487 Analysis 


paar wos 


The WE cutput froa tne SSFFT model iS i ic 
units, waick is a probiem when coupariny differext todeis. 
However, for this investigation the proplem was mininal. As 
expressed by equation 2.6, the reference wavenunper for 100 
a2 with a reference sound speed of 4921.26 ft/sec (1509 
m/sec) is 0.1277 1/ft. After correction of the reference 
wavenumber for a range increment of 52.3 ct and the 2045 
points used to represent the wavenunber spectrun, the 
Maximum beta value should be (@.1243 1/ft which agrees 
reasonably well with the maxigua beta vaiues in figures 3.3, 
3.4, and 3.5. Contrary to theory the nuils ar2 not equally 
spaced and the expected delta beta value of 9.001592 1yit 
(equation 2.19) is net represented by any of tae null spac- 
ings. Stamey [Ref. 3], observed cae saze non-uniftora 
Spacing of beta nulls. The reasoa way the source depth is 
not represented by any of the nuli spacinys may be traced to 
tae improper aull placement in the raage domaix. Tae Lack 
of eguai spacing between tie nuils nust be attributed to 
something elise, possibly the pressure wave was inproperliy 
represeated as a continuous signal or nuagerical errors exist 
within the source ccde. Whatever the reason for this 
gifference, it must Le resolved before the SSFPT can be usei 
with tne WI to deterzine the source iepth. 
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C. IFD 


The IFD model [Ref. 12] wWiicnh was iastalicu aac 
tested in September 1983 at NPS [Ref. 13] was the next cecel 
chosen aS a sound source for the ki. Tais decision was rads: 
kecause the IFD was capable or properiy handling the ocean 
bottom interactions and was the only otaoer acoustic PE nocel 
an residence at NPS. A copy of the source code Listing is 
provided in appendix ft. 

The IFD also presented fprodieas which will se 
discussed further. The same initial test ran used witn the 
SSFFT waS used with the IFD. Since the IFD requires rore 
bottom infornation, additional testin5 was perftorzed in 
order to obtain a fully absorbing potton. Fijure 3.6 iilus- 
trates the environment as seen py the IF) nodcei. 

The densities in the water mass, sedinent, and arti- 
ficial attenuation layer were set to a constant of 1.0 and 
the attenuation in the sedjinent was adjusted until the ¢iia- 
ination of bottom interference was observed. Figure 3.7 is 
an example of the IFD transaissioa loss output for a 500 
source/receiver depth, a water depta of 2000 a, aa uprer 
surface of the artificial attenuation layer at 3000 a, ania 
dower pressure release surface at 400) a. The fottoz atten- 
uation used to produce figure 3.7 was J.0916 a3/wavelength. 
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Figure 3.6 IFD Environment. 
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The location 9f the nulls in figure 3.7 coiuciie witn t 

predicted by equaticn 2.4 a1f the higa rrejvency noise is 
filtered out by either a running averaje process or low sés: 
digital filter. A careful examination of the transadission 
loss curve (figure 3.8) reveals the real problen, which is 
the long wavelength noise near tae wavelength of i. terest. 
The long wavelength noise will cause erroneous spikes on tie 
WI plots after being Fourier transformed. Coatinued anal- 
ysis revealed that the adise could be eliminated if tke 
thicknesses of the artificial attenuation layer ani the 
bottom pressure release surface are caosen so that they are 
at least twice the water depth. Ine choice of a snail 
vertical grid step will reduce the thickness reguired for 
the pressure reiease surface, but this version of the IF) 
progran is linited to 5000 vertical grid points due to soit- 
ware restrictions. Figure 3.9 is a transmission loss curve 
obtained from a 500 m source/receiver depth, a 20900 a water 
depth, a400) om upper level artificial attenuation layer, 
and an 3000 m bottom pressure release surface. Comparison 
of figures 3.1 and 3.9 wits the classical Lioyd's mirror 
transmission loss (figure 2.2) snows that tne shape of the 
curve beyond 30 kn is greatiy improved with the IFD. The 
number of vertical grid steps was 5090. when the IfFD is 
used in a shallew water situation and at 25 hz, tne noise 
does not appear td be present. Figure 3.10 is a plot of the 
transmission loss for a signal at 25 nz, witn the source/ 
receiver depta at 50 m, the bottom sloping up froa 350 to 50 
DB, and an artificial attenuation iayer and bottom pressure 
release surface located at 750 and 1000 a, raspectively. 
Toe noise whith was so evident in rigure 3.7 is absent froa 
the Shallow sater case. Producing tae noise free curves 
requires excesSive computer time except for shallow water 


cases. 
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To reduce the amount of the oressure infornatcicn ana 


es 


yet preserve the detaiis, the IFrDdD tests were contucte: 
between 10 and 40 km in range. The reference wavenudber for 
tne test input with the IFD was 0.4139 1/a wWhicr, aiter 
correcting for a cange increment of 7.5 a witn 4)9€ peints, 
produced a bveta maximum of 0.4188 1/n. Figures 3.11, 3.12, 
and 3.13 indicate agreenent with the beta maxinum ana a 
delta beta of 0.00623 1/m can be found agony tie Ruii Spacr- 
ings. when the vertical grid step is reaucead in size, tne 
output signal contains amore iaforzation of finer détail. 
The finer detailed irformation will, iz noise is preseLt, 
produce more clutter on the WT depiction (figure 3.13). 5 
Numerous WI depictions a peak was observed at the ainiauc 
and maximum waveauanbers. Stamey celieved that tre rignt and 
left intensity maxina corresponded t9 bean €levation anjics 

£ 0 and 30 degrees, respectively. de further stated that 
the lert naxina aiso may be relatei to aigsorithnic inaccira- 
cies [{Ref. 3}. The “U-shaped phenoazaon" cbhserved auria; 
Staaey's investigaticn or the SSirft, is present witt tae IF) 
and the nulls are not equally spaced. 
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~aengaenon™ appears to cc 


Tne reason f9r the wyespayed 
related to noise in tae envelope function yield wnic. 25 
further appliried ahen tke Hankel function ts insertsi >? 
obtain tne conplex pressure Using a Low pass fiit epee 
jner easing the vertical grid step size appear to reduce the 
problen. In figure 3-8, the nulis occir at 33.30 ane 33.3% 
km when 4 single null should exist at 33.33 kB. This condi~ 
tion or noise produces the wost JZevastatiag effects 57 
maskin, the desired gGensity spectruda. A narrow band pass 
fiiter would probably solve this problen in the short tera 
but eventualiy the IFD will require Earther study to eliou= 
nate the noise. dhen the shallow case is used is the at 
(figure 3.174), the ny-sha ped phenomenon" is acsent ana tne 
null locations are more distinct. Howe veLls in the spaliow 
water case, the ootton interaction produces nigh Ere quescy 


interference but tne sulil spacings occur near tne expected 


walue ror the source depth (0.0628 1/) - 
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Two fparaboiic eguation models, tae SSFFT and IiD, were 


used to predict sound fielis for comparison with the Lioyd's 
micror interference pattern in the range domain and tnen in 
the wavenumber domain. In the range domain, tie SSFFT 
“aproperiy handled the location of pressure nuils ané 
Gisplayed bottom interrerence. Wnile these weaxnesses in 
the SSFrT were absent in the IFD, the inadvertent and dele- 
terious insertion of noise was extensive. In tae wavenunter 
domain, tne noise in the IFD was prohibitive in ail but the 
deep sediment and shaliow water cases. Tha IFD producei 
results which were ccmparable with the SSFFI when the noise 
was not a factor. The range step and tae bottom attenuation 
algoritnms of the SSFFT should be reviewed in orier to 
correct those problems descriked avsove. aa the IF), the 
noise prodlem will reguire further investigation. fs tac 
problem in tne IFD is not corrected, this nodel will se 
severely liaited. With the guidance provided by Lauer 
(Ref. 2], the wavenumber domain can provide tae detailed 
information needed to quickly correct inconsistencies 
opserved between acoustic models. Therefore, the WT should 
be investigated further to ensure taat any soitware weak-~ 
nesses are eliminated and a thorough understacding, of the 
process is assured. 

Experiments in the operational environment should be 
plaoned to test the WI with actual sound sources. The ship- 
board experiment could be easily accomplished since the WT 
can ce implemented on micro-computers if the pressure infor- 
Mation is provided from an outside source. The wi, if it is 
going to ke considered asa metnaoad to determine source 
depth, appears to have three major saortcomings: 
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1. A clean C# signal with phase iaforaation is ese 22- 


tial, which will require guaarature deaocdalation one 
extensive filtering. 

2. Prediction models are required which will proviie 
egualiy spaced beta nulls so taat software can £ce 
developed to automate the process. 

3. High resolution processing eguipment is necessary 
because the change in beta spacing for shallow sound 


sources is small (figure 2.5). 
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